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ABSTRACT. - A preliminary morphometric study based on 17 landmarks was carried out, using 104 specimens of six spe¬ 
cies of sea catfishes from both South American coasts, the four accepted species of the genus Bagre and two reference spe¬ 
cies ( Ariopsis seemanni and Arius platypogon). The landmarks were mainly located on the head and at the origin and end 
of fins. The partial warps weighted matrix was used to perform the discriminant analysis and the species pair-wise compari¬ 
sons in order to establish if there were differences between the species. Thereafter, a partial-warps weighted matrix was 
generated with the average species consensus. Using PHYLIP maximum likelihood routine (CONTML), an additive tree 
was obtained. Shape changes were described over the tree using TPS tree (square change parsimony). Bagre marinus 
(western Atlantic) seems to be the most basal species of the genus. The species showing the larger shape changes are B. 
pinnimaculatus (eastern Pacific) and B. bagre (western Atlantic). Those two species are apparently more closely related to 
each other than to the other two species of Bagre. 


RESUME. - Etude systematique du genre Bagre (Siluriformes : Ariidae). Approche morphometrique. 

Une etude morphometrique preliminaire, basee sur 17 reperes, est realisee sur 104 specimens appartenant a six especes 
de poissons-chats marins des deux cotes sud-americaines : les quatre especes acceptees du genre Bagre et deux especes de 
reference ( Ariopsis seemanni et Arius platypogon). Les reperes sont principalement localises sur la tete ainsi qu’a l’origine 
et a l’extremite posterieure des nageoires. Afin de mettre en evidence les differences interspecifiques, une analyse discrimi- 
nante et des comparaisons interspecifiques deux a deux ont ete realisees sur la matrice des flexions partielles. Ensuite, une 
matrice des flexions partielles a ete construite a partir du consensus obtenu sur les poissons moyens par espece. Un arbre 
par addition a ete obtenu en utilisant la routine du maximum de vraisemblance du progiciel PHYLIP (CONTML). Les 
changements de conformation ont ete decrits d’apres l’arbre par utilisation de TPSTree (parcimonie sur les moindres car- 
res). Bagre marinus (Atlantique ouest) semble etre l’espece la plus plesiomorphe du genre. Les especes montrant les plus 
importants changements morphologiques sont B. pinnimaculatus (Pacifique est) et B. bagre (Atlantique ouest). Ces deux 
especes sont apparemment plus etroitement apparentees Tune a Tautre que les deux autres especes de Bagre. 
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The family Ariidae (order Siluriformes), commonly 
known as sea catfishes, is widely distributed around the world 
in all the tropical and subtropical continental shelves, but the 
intrafamilial phylogenetic relationships are poorly known 
(Kailola and Bussing, 1995; Acero, 2002). Ariidae includes 
between 130 (Acero and Betancur, 2002) and 200 species 
(Acero, 2003) that inhabit mainly in estuaries and turbid 
coastal waters, but some species occur below 100 m depth 
and others dwell basically in freshwater. However, the num¬ 
ber of valid genera is unclear as well as their distribution. 
Betancur et al. (2003) made a preliminary analysis of the sys¬ 
tematic of the family in both coasts of northern South Ameri¬ 
ca, taking into account some osteologic and external charac¬ 
ters. The resulting phylogenetic hypothesis is a cladogram 
where they show that there are many inconsistencies in the 
generic and suprageneric systematic of the family. 


Simple morphometric procedures have been established 
and applied in systematic and taxonomic studies for more 
than a century, being employed with greater frequency since 
the decade of 1970, due to the development of more elaborat¬ 
ed methods (Chapman et al., 1997). Shape analysis, morpho¬ 
metries, and statistics are available tools for researchers in 
order to increase the foundations of their results (Chapman, 
1990). Morphometric techniques have been used to establish 
limits between species, populations, and stocks (Stracahan 
and Kell, 1995; Cadrin and Friedland, 1999; Littman et al., 
2000), as well as to describe biogeographic patterns and infer 
the phylogeny of a group (Loy and Corti, 1996). 

Rohlf (1990) defined morphometries as the quantitative 
description, analysis and interpretation of shapes and their 
variation in biology. He determined that description and 
comparison of structures’ shapes are necessary in any sys- 
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tematic study that is based on organisms morphology. Sev¬ 
eral studies have shown that the recently developed tech¬ 
niques of geometric morphometries are extremely powerful 
descriptive tools (Swiderski et al., 2000). 

However, the use of the morphometric techniques is still 
under debate. Swiderski et al. (1998,2000) and Zelditch et al. 
(1998) consider that the homology principle may be applied to 
the relative deformation warps, which are used as characters 
to infer phylogenetic hypothesis. Adams and Rosenberg 
(1998), Rohlf (1998b), Monteiro (2000), and Roth and Mercer 
(2000), on the other hand, opine that the mathematical proper¬ 
ties of the warps lack biological meaning; therefore, the use of 
those data as cladistic characters may lead to erroneous inter¬ 
pretations in the reconstruction of the evolutionary history of 
a group. Nevertheless, the group of warps coming from a 
selection of anatomical landmarks may be useful to infer a 
phylogenetic, non cladistic, hypothesis based on a technique 
such as square change parsimony and maximum likelihood 
because those procedures are invariant to the rotations inher¬ 
ent to the shape space (Rohlf, 2002). 

The purpose of this study is to demonstrate how shape 
has changed from one species to another and that morpho¬ 
metric techniques could be used to infer phylogenetic rela¬ 
tionships as well as congruent hypothesis of classification. 
The systematics of the genus Bagre, a sea catfish group 
endemic to the Neotropics and represented by four species, 
two from the Caribbean (B. bagre and B. marinus ) and two 
from the eastern Pacific (B. panamensis and B. pinnimacula- 
tus ), are herewith analyzed. 


MATERIAL AND METHODS 

A total of 104 specimens representing six species were 
studied. Additionally to the four species of Bagre, two spe¬ 
cies from the eastern Pacific were chosen as the ‘outgroups’ 
for shape analysis. They were Ariusplatypogon, considered 
by Betancur et al. (2003) as the sister species to the genus 
Bagre as a whole, and Ariopsis seemanni, selected because it 
is included in the same suprageneric group as Bagre (Betan¬ 
cur et al., 2003). The Caribbean specimens were collected in 
the Guajira peninsula (Colombia) (11° 25’N, 73° 5’W) and 
in Margarita Island (Venezuela) (11°N, 64°0rW); the Pacif¬ 
ic species were all collected in Buenaventura, Colombia (3° 
55’N, 77° 05’W). 

The specimens were preserved in formaline for two to 
five days, in order to avoid bending and possible errors 
caused by differences in the position of selected landmarks. 
The digital images were obtained using a flat bed scanner 
(Epson ESI200 C). The images, taken only on the right side 
of the fishes, were carefully taken to not mask important 



Figure 1. - Location of the 17 selected landmarks (black dots) over 
one image of Bagre marinus species. [Localisation des 17 reperes 
selectionnes (points noirs) sur Bagre marinusj 


structures, and were edited in the program Adobe Photoshop 
v.3.0. Finally, 17 landmarks were digitalized employing the 
TPSdig v.1.39 application developed by Rohlf (2003) 
(Fig. 1). 

The shape variation was analyzed using TPSrelw vl.23 
application (Rohlf, 1998c). The technique consists in the 
superimposition and alignment of all the specimens to a 
consensus, showing the changes in function of the aver¬ 
aged Procrustes distance in Kendall’s shape space among 
the consensus and the entire sample of data, using general¬ 
ized orthogonal least-squares (GLS) (Rohlf and Slice, 
1990). The alignment is then projected to a tangent shape 
space (Rohlf 1999b, 2000) in which a principal compo¬ 
nents arrangement known as relative warps is performed 
using a partial warp matrix (Bookstein, 1989, 1991, 1996; 
Rohlf, 1996, 2000). Relative warps are just principal com¬ 
ponents vectors in the space used to describe the major 
trends of the shape variations between the specimens as 
deformations to a reference (Rohlf, 1999a). TPSsmall 
{ibid) application was used to check if the tangent projec¬ 
tion is close to the original Kendall’s shape space defined 
by the whole data. 

The weighted matrix of the partial warps that is obtained 
in TPSrelw was introduced in the statistic software SYSTAT 
v.9.0 to perform a step by step discriminant analysis using 
jackknife procedure. Once it was establish that there were 
significant statistic differences among the species, a centroid 
image from all the species was generated with TPSsuper 
(Rohlf, 1998a), as a value to estimate a matrix of the species 
and generate a weighted matrix of the partial warps to the 
averaged value of the landmarks configuration for each spe¬ 
cies with TPSrelw (Rohlf, 2000). 

Finally, a tree was estimated using the maximum likeli¬ 
hood procedure developed by Felsenstein (2001) in the pro¬ 
gram PHYLIP v 3.6 using the routine CONT (continuous 
variables, no rotted tree, randomized order of the sample, 
and 100 replicates). This tree produces a phylogenetic 
hypothesis. After that the tree is introduced in the program 
TPStree (Rohlf, 1999a) to find how its topology can reflect 
the changes in form. 
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Table I. - Maximum, minimum and medium values of the Procrust¬ 
es and tangent distances. [Valeurs maximales, minimales et moyen- 
nes des distances Procustes et des distances tangentes.] 



Procustes 

distances 

Tangent 

distances 

Minimum 

Maximum 

Medium 

0.022922 

0.186963 

0.089733 

0.022908 

0.186690 

0.089690 


Table II. - Statisticals for the discriminant analysis. [Statistiques 
pour l’analyse discriminate.] 


Statistical 

Value 

F Value 

Probability 

Wilkis’s lambda 

0 

42.140 

0 

Pillai’s trace 

4.19 

24.436 

0 



Figure 2. - Tangent against Procrustes distances regression for all 
analyzed specimens. [Regression entre distances dans I’espace 
tangent et distances Procustes pour tons les specimens analyses.] 

RESULTS AND DISCUSSION 

Tangent distances vs. Procrustes distances 

Procrustes distance and tangent space distance of the 
shapes of the sea catfishes were compared. The relation of 
the regression between the shape spaces was almost lineal, 
with none of the studied specimens showing any significant 
deviation from the slope determining the relation between 
the shape spaces (Fig. 2). Table I includes the maximum, 
minimum, and mean values of the Procrustes and tangent 
distances, as well as the slope (0.999439) and correlation (1) 
values. Marcus et al. (2000), in their morphometric analysis 
of mammal skulls, show slope and correlation values of 
0.9947 and 0.9997, respectively. Even though there are dif¬ 
ferences in the configuration of the spaces, the distortion of 
the multidimensional space projection (Procrustes distance) 
in the tangent space is small enough to allow the use of the 
traditional statistic methods to describe the variation in 
shape, as it was proposed by Rohlf (1999b). 



Figure 3. - PCA Ordination plot showing the 104 specimens over 
two first relative warps. A: Relative warp lx vs ly for 2D data; 
B: relative warp lx vs 2x for 2D data. [Analyse en composantes 
principals montrant les 104 specimens pour les deux premieres fle¬ 
xions relatives. A : Flexion relative lx enfonction de ly ; B : Fle¬ 
xion relative lx enfonction de 2x.[ 

Relative deformation warps and discriminant analysis 

After aligning and overlapping the group of landmarks 
for the 104 specimens, the first three deformation axes accu¬ 
mulate 80.5% of the information on the body form variation. 
The graphics of the relative deformation warps show that the 
species get organized in three groups in the ordination 
(Fig. 3A: lx vs ly; Fig. 3B: lx vs 2x), which are relatively 
well separated. Bagre bagre (points 35 to 51 in the ordina¬ 
tion) and B. pinnimaculatus (points 88 to 104) appear in the 
positive space of lx warp. Bagre marinus (points 52 to 65) 
and B. panamensis (points 66 to 87) appear in the middle 
part of lx warp. On the other hand, Ariopsis seemanni 
(points 7 to 34) and A. platypogon (points 1 to 6) appear in 
the negative space of lx warp. Finally, B. marinus and B. 
panamensis as well as A. platypogon and A. seemanni are 
clearly separated on the 2x warp (Fig. 3B). 

The results obtained with the discriminant analysis lead to 
know the existence of statistical differences between the stud¬ 
ied species (Tab. II). With the Mahalanobis distance classifi¬ 
cation, 100% of the individuals are correctly assigned to their 
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Table III. - Classification matrix after Jackknife procedure in the 
discriminant SYSTAT routine. Rows correspond to species and col¬ 
umns correspond to specimens correctly assigned (%) to their origi¬ 
nal taxa. [Matrice de classification apres procede de Jackknife dans 
la routine discriminante SYSTAT. Les lignes correspondent a 
l ’espece et les colonnes correspondent aux specimens correctement 
assignes (%) a leurs taxons originaux.] 


Species 

A. platypogon 

A. seemanni 

B. bagre 

B. marinus 

B. panamensis 

B. pinnimaculatus 

% assignation 

A. platypogon 

6 

0 

0 

0 

0 

0 

100 

A. seemanni 

0 

28 

0 

0 

0 

0 

100 

B. bagre 

0 

0 

17 

0 

0 

0 

100 

B. marinus 

0 

0 

0 

14 

0 

0 

100 

B. panamensis 

0 

0 

0 

0 

22 

0 

100 

B. pinnimaculatus 

0 

0 

1 

0 

0 

16 

94 

Total 

6 

28 

18 

14 

22 

16 

99 


relative species; after the jackknife procedure the percentage 
diminish to 99% (Tab. III). With the canonic functions, the 
cumulative variation for the first and second factor is 85.66%. 
The ordination of the species using the first three canonic vari¬ 
ables was similar to the one shown in figure 3; that is each one 
of the species pairs (A. platypogon-A. seemanni, B. marinus- 
B. panamensis, and B. pinnimaculatus-B. bagre ) shows a 
shape relatively similar within the space (Fig. 4). 

The way that the species are placed in the ordinations 
(Figs 3,4) suggests that there is a direction in the change of 
the body shape that may be interpreted as a signal of the phy- 
logeny of the group, due to the low diversity of shapes in the 
bulk of analyzed data. This differs clearly to what was found 
by Marcus et al. (2000) in their mammal skulls study. 
Zelditch et al. (1995), Marcus et al. (2000) and Swiderski et 
al. (2000) have clarified the objections about the homology 
principle (de Pinna, 1991) with this 
kind of data sets. Therefore, the 
features analyzed may be used to 
infer phylogenies without using 

*4 

Figure 5. - Schematic shape changes 
over first relative warp. A: Positive 
lx side; B: Negative lx side. Dots 
shows consensus configuration, arrows 
shows direction and magnitude of the 
vectors of the shape change. [Change- 
ments schematiques de conformation 
enfonction de la premiere flexion rela¬ 
tive. A : Coordonnees lx positives ; 

B : Coordonnees lx negatives. Les 
points representent la conformation 
consensus, lesfleches representent la 
direction et la norme des vecteurs de 
deformation.] 



Figure 4. Plot of the first two canonical factors, after Jackknife 
SYSTAT routine, showing specimens clusters, using species as 
dummy group variant: ( i ) Bagre pinnimaculatus , (+) B. bagre , 
() B. marinus (▼) B. panamensis , (x) Ariopsis seemanni , and 
( •) Occidentariusplatypogon. [Deuxpremiers facteurs canoni- 
ques, d'apres la routine Jackknife SYSTAT, montrant les regroupe- 
ments de specimens, en utilisant les especes comme groupes d’ob- 
jets dijferents.] 

cladistic methods because the codification of continuous 
variables is unjustified (Felsenstein, 1988). 

Ariusplatypogon and A. seemanni are located in the ordi¬ 
nations in a position opposite to the species of Bagre, which 
may be in agreement with the outgroup principle. Those taxa 
(i.e., species, genera) that may be related but are not included 
in the target taxa may be considered part of the outgroup. 
Ideally the outgroup should be the sister group to the ingroup 
(Brooks et al., 1984), which then may give an indication of 
the shape changes and infer evolutionary novelties. The spe¬ 
cies closer to the consensus are B. marinus and B. panamen- 
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Figure 6. - Schematic shape changes 
over second relative warp. A: Positive 
IX side; B: Negative IX side. Dots 
show consensus configuration, arrows 
show direction and magnitude of the 
vectors of the shape change. [Chan- 
gements schematiques de conformation 
en fonction de la deuxieme flexion 
relative. A : Coordonnees lx positives ; 
B : Coordonnees lx negatives. Les 
points representent la conformation 
consensus, les fleches representent la 
direction et la norme des vecteurs de 
deformation.] 



Figure 7. - Phylogenetic hypothesis tree (additive tree), using 
unwarped average species from TPSsuper procedures, as hypothetic 
taxonomic units for each species. Weighted matrix (TPSrelw) was 
generated with average specimens; tree was generated using Phylip 
CONT routine (maximum likelihood, random order samples, no rot¬ 
ted tree, 100 replicates). [Hypothese phylogenetique (arbre additif), 
a partir des consensus moyens des especes (TPSsuper) utilises 
comme representants hypothetiques pour chaque espece. La matrice 
des flexions partielles a ete creee a partir du consensus des speci¬ 
mens et l’arbre est calcule avec la routine CONT de PHYLIP (maxi¬ 
mum de vraisemblance, ordre aleatoire des echantillons, arbre non 
enracine, 100 replicats).] 

sis , while B. pinnimaculatus and B. bagre are placed as the 
most distant to the consensus (Figs 3,4). 

Description of the changes in shape in relative 
deformation warps 

The most evident changes in shape on the first warp 
appear in: 

A. Head, mainly in the region anterior to the snout (land¬ 
marks one to three) (Fig. 5). Landmark 1 may elongate, 
causing the forward displacement of the insertion of maxil¬ 
lary barbel (landmark 2) and the anterior border of eye (land¬ 
mark 3). The movement of landmark 3 corresponds to the 
increase in eye diameter. 

B. Dorsal fin (landmarks six and seven); dorsal fin base 
length remains unchanged, however predorsal fin length 
changes (Fig. 5). 


C. Anal fin (landmarks 13 and 14); the length of the fin 
base changes, however its end remains invariant. Due the 
movement of the anterior extreme, the fin moves forward and, 
as a consequence, the pelvic fin position changes as well. 

D. Insertion and orientation of pectoral fin (landmarks 16 
and 17); the insertion angle varies causing that the complex 
moves ventrally to end up parallel to the longitudinal body 
axis (Fig. 5). 

In the second relative deformation warp the body shape 
changes are associated to a dorso-ventral displacement of the 
landmarks. Therefore, when the shape space gets depressed, 
anal fin base length and adipose fin base length increase. If 
the body gets deeper, the base lengths of the two fins decrease 
(Fig. 6A, B). 

Additive trees and shape changes 

Figure 7 shows the additive tree obtained using PHYLIP, 
on which the variation of the landmarks position in the shape 
space on each one of the tree nodes was studied. When the 
distance scale is analyzed, it can be seen that B. pinnimacu¬ 
latus and B. bagre are widely separated from the other two 
species of the genus. Furthermore, it is clear that the land¬ 
mark configuration in B. marinus is the most plesiomorphic 
of the genus, making it the most adequate for shape compari¬ 
sons. 

Phylogenetic hypothesis 

The final result, the tree showing the relationships exist¬ 
ing between the species of Bagre and the outgroup species, 
is presented in figure 7. Basally the tree is a trichotomy, with 
the relationships between A. platypogon, A. seemanni , and 
Bagre unresolved. Therefore, the evidence provided by the 
morphometric characters was not sufficient to determine 
which of the two outgroup configurations is closer to Bagre. 
The tree shows the relationships between the species of 
Bagre in the same way as the ordinations (Fig. 7). 
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Table IV. - Anal fin rays range in the genera and selected species of 
neotropical sea catfishes. It is indicated whether the taxon occurs in 
the western Atlantic (A) or in the eastern Pacific (P) and the number 
of species in each ocean. Cathorops n.sp. is in description process 
(Betancur and Acero, in press). [Nombre de rayons de la nageoire 
anale chez les genres et les especes selectionnes de poissons-chats 
des mers tropicales. II est indique si le taxon apparait dans VAtlan- 
tique occidental (A) ou dans le Pacifique oriental (P) et le nombre 
d’especespour chaque ocean. Cathorops n. sp. est en cours de des¬ 
cription (Betancur et Acero, souspresse).] 


Genera or species 

Anal fin rays 

A 

P 

Ariopsis 

17-20 

3 

2 

Arius 

17-23 

6 

4 

A. lentiginosus 

23-28 

- 

1 

A. platypogon 

18-20 

- 

1 

Bagre bagre 

29-37 

1 

- 

B. marinus 

22-28 

1 

- 

B. panamensis 

25-30 

- 

1 

B. pinnimaculatus 

27-32 

- 

1 

Cathorops 

19-25 

9 

4 

C. multiradiatus 

24-28 

- 

1 

Cathorops n.sp. 

21-25 

1 

- 

Galeichthys 

15-18 

- 

1 

Selenaspis 

17-21 

3 

1 


For the phylogeny of the genus Bagre, all the available 
recurrent evidence produced by geometric morphometry 
(relative deformation warps, discriminant analysis, and addi¬ 
tive trees) points towards the strong similarity between B. 
bagre and B. pinnimaculatus. They are also placed as the 
most derived species within the genus, which tends to back 
the systematic hypothesis of their close relationship. On the 
other hand, it seems that B. marinus and B.panamensis are 
the most generalized species, with B. marinus being the most 
plesiomorphic and, therefore, hypothetically the oldest liv¬ 
ing species of the genus. 

Another feature that seems to support the proposed 
hypothesis is the evolution of the number of anal-fin ele¬ 
ments. If our hypothesis is correct, the number of anal-fin 
elements has been increasing in the same direction than the 
evolution within the genus Bagre. Hence, B. marinus, con¬ 
sidered the most basal species, has 22-28 anal-fin rays, B. 
panamensis has 25-30, B.pinnimaculatus has 27-32, and B. 
bagre has 29-37. The other species of sea catfishes in the 
western Atlantic and the eastern Pacific have 15-25 anal-fin 
elements (except Arius lentiginosus, 23-28, and two 
Cathorops species, 24-28) (Cervigon, 1991; Kailola and 
Bussing, 1995; Marceniuk, 1997; Acero, 2002; Acero and 
Betancur, 2002; Betancur and Acero, in press) (Tab. IV). 
Therefore, it seems that the speciation processes linked to 
the evolution of Bagre have been accompanied by an 
increase in the number of anal-fin elements. The species of 


Bagre have more anal-fin rays than any catfish of the clade 
Ariopsis - A. platypogon - Bagre - Selenaspis, which has 
been considered monophyletic (Betancur et al., in press). 
Since no species of this clade, other than the members of 
Bagre, has more than 21 rays, it seems that the origin and 
evolution of the genus has been linked to an increase in the 
number of anal-fin elements. 
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